Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 September 2011 (MN 1^1^ style file v2.2) 



Anisotropies in the Diffuse Gamma-Ray Background from 
Dark Matter with Fermi LAT: a closer look 



o 

Oh' 



^' 
o 

(N 

oo 

p 

o 
o 



X 



A. Cuoco^'^, A. Sellerholm^, J. Conrad^ and S. Hannestad^ 

^ The Oskar Klein Centre for Cosmo Particle Physics, AlbaNova, SE-106 91 Stockholm, Sweden 

^Departm,ent of Physics and Astronomy, University of Aarhus, Ny Munkegade, Bygn. 1520 8000 Aarhus Denmark 



7 September 2011 



ABSTRACT 

Wc perform a detailed study of the sensitivity to the anisotropies related to Dark 
Matter (DM) annihilation in the Isotropic Gamma- Ray Backgromid (IGRB) as mea- 
sured by the Fermi Large Area Telescope (Fermi-LAT). For the first time, we take into 
account the effects of the Galactic foregrounds and use a realistic representation of 
the Fermi-LAT. We implement an analysis pipeline which simulates Fermi-LAT data 
sets starting from model maps of the Galactic foregrounds, the Fermi resolved point 
sources, the extra-galactic diffuse emission and the signal from DM annihilation. The 
effects of the detector are taken into account by convolving the model maps with the 
Fermi-LAT instrumental response. We then use the angular power spectrum to char- 
acterize the anisotropy properties of the simulated data and to study the sensitivity 
to DM. We consider DM anisotropies of extra-galactic origin and of Galactic origin 
(which can be generated through annihilation in the Milky Way sub-structures) as 
opposed to a background of anisotropies generated by sources of astrophysical origin, 
blazars for example. We find that with statistics from 5 years of observation Fermi is 
sensitive to a DM contribution at the level of 1%-10% of the measured IGRB depend- 
ing on the DM mass TOj^- and annihilation mode. In terms of the thermally averaged 
cross section (uj^v), this corresponds to ^ 10~^^cm'^s~-^, i.e. slightly above the typi- 
cal expectations for a thermal relic, for low values of the DM mass m^ < 100 GeV. 
The anisotropy method for DM searches has a sensitivity comparable to the usual 
methods based only on the energy spectrum and thus constitutes an independent and 
complementary piece of information in the DM puzzle. 



1 INTRODUCTION 



Combined studies of the cosmic microwave background ra- 
, diation, supernova cosmology and large Galaxy redshift sur- 
veys provide nowadays compelling evidence of the existence 
of non-baryonic Dark Matter (DM) (|Komatsu at al. 20081 

|Tegmark et al. 2006[ ). Consistency with the observed struc- 
ture of the Universe, especially at large (above galactic) 
scales favors Cold Dark Matter (CDM), i.e. the particles con- 
stituting the cosmologically required dark matter have to be 
moving non-relativistically. This property is always fulfilled 
by particles of mass in the GeV to TeV region that interact 
with the weak interaction strength, i.e. weakly interacting 
massive particles (WIMPs). WIMPs can annihilate or decay 
to detectable standard model particles, in particular gamma- 
rays, giving rise to the possibility of "indirect detection" 
(|Bertone et al. 20041 [Jungman et al. 1995||Bergstrom 2000[ 

IBaltz et al. 2008|) . A standard approach is to look for spec- 
tral signatures in the energy spectrum of gamma-rays from 
celestial objects, either targeting the continuum spectrum 
produced by decay of secondary pions or the more feeble 
signature of a line due to DM annihilating directly into pho- 
tons. 

Another possible experimental signature is given by 



differences in the spatial distribution of gamma-ray signals 
induced by DM as compared to conventional astrophysical 
sources. In particular, the small scale fluctuations in an 
otherwise isotropic gamma-ray background (IGRB) might 
be different in the case of annihilating DM, since its flux 
scales as the squared density. After the pioneering study 
of Ando and Komatsu (|Ando fc Komatsu 2006|) , followed 
by two other more detailed works (|Ando et al. 20061 
I Ando et al. 2007)) . the issue of DM annihilation and 
anisotropies in the isotropic gamma background has raised 
increasing interest. In connection to the extra-galactic 
gamma-ray background, anisotropy has been further stud- 
ied in (|Cuoco et al. 20061 ICuoco et al. 20071 ICuoco 2008)l 
while predictions from the Millenium II simulation have 
recently been described in Zavala et al. (2009). It has also 
been realized that Galactic substructures can produce an 
almost isotropic gamma background with similar fluctuation 
properties as the extragalactic background and with promis- 
ing chances of detection ([Siegal-Gaskins 2008[ [Ando 20091 



[Siegal-Gaskins and Pavlidou 2009[ [Hensley et al. 2009|); 



see also (|Pieri et al.lOOTi IPieri et al. 2009|) ). Other 
works have investigated the anisotropy pattern result- 
ing from both the Galactic and extragalactic contri- 
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bution (|Fornasa et al. 20091 [Hooper fc Serpico 20"07l 
llbarra et al. 2010)) or from the the DM annihilation 
around intermediate mass black holes (|Taoso et al. 200"9)) . 
Anisotropics from DM annihilation appearing in the radio 
sky have been investigated in ( [Zhang and Sigl 2008[ ). 
Finally,- apart from the power spectrum - an interesting 
approach is to investigate the probability distribution of 
the fluctuations in order to detect possible departures 
from Gaussianity (or, rather, from Poisson statistics, in 
the common situation of low number of counts), thereby 
providing another handle on the separation of conven- 
tional astrophysical and DM contribution (|Lee et al. 20081 
IDodelson et al. 2009|) . Anisotropy properties of astro- 
physical processes other than DM contributing to the 
isotropic gamma background include resolved point 
sources (!Ando c t al. 2006)) . galaxy clusters and shocks 
l[Miniati e t al. 2 0071 IKeshet et al. 2003)) and emission from 
normal galaxies (|Ando fc Pavlidou 2009|) . 

The Fermi Large Area Telescope (Fermi-LAT), which 
was launched June 11, 2008, is the most suitable tool 
for these kind of studies due to its large field of view 
(~ 2.4 sr), energy coverage from 100 MeV to 300 GeV 
l[Atwood et al. 2009)1 and very good angular resolution. We 
will thus focus in the following on the prospect to detect 
DM anisotropies with Fermi-LAT. 

All of the above studies assume either that the impact 
of Galactic foregrounds on the DM induced anisotropies is 
negligible or that the Galactic foregrounds can be cleaned at 
a suitable level to allow DM anisotropies studies. Further, 
the effects of the Fermi-LAT detector are generally taken 
into account only in an approximate way, though a detailed 
representation of the Fermi instrument is publicly available. 

The aim of this work is to address the above 
issues in detail. In particular our approach will 
be the following: we include the Galactic fore- 
grounds explicitly using a model based on the 
GALPROP code ( [Strong and Moskalenko 1998[ 

[Moskalenko and Strong 1997| [Strong et al. 2004b[ ) aiid 
we then use this model to build masks to select the most 
promising regions in the sky suitable for the analysis of the 
anisotropies. In these regions we study the anisotropies in 
the scenario in which a DM contribution is present over the 
Galactic foregrounds and the extra-galactic emission from 
blazars compared to the case without DM and we derive 
the sensitivity to the DM signal employing anisotropy 
studies with Fermi-LAT data. Both the cases of Galactic 
DM in substructures and extra-galactic DM are considered 
separately. The effect of the Fermi-LAT response, including 
the effects of the point spread function, the effective area, 
as well as non-uniformities in exposure, are included in 
all calculations at a realistic level by using the instrument 
response functions provided by the Fermi Science Tools. 

This approach is suitable for a sensitivity study like the 
one presented here, and does not aim to fully reproduce a 
real data analysis for which foreground cleaning (not per- 
formed in this study) most likely will be required. One of 
the findings of our analysis is indeed that even with a very 
generous masking, foregrounds still have a sizeable contribu- 
tion to the overall anisotropies and represent an important 
limitation to DM sensitivity. This implies that foreground 
cleaning will be a crucial step in the data analysis pipeline. 
On the other hand, if cleaning can be achieved in an effec- 



tive way, this also means that the sensitivity to DM can be 
further improved with respect to the results derived in this 
paper. Our results can thus be regarded as conservative from 
this point of view. Foreground cleaning will be the subject 
of forthcoming investigations. 

The paper is structured as follows. In section II we in- 
troduce and discuss the various components contributing to 
the Galactic and extragalactic gamma-ray emission. In sec- 
tions III and IV we review their anisotropy properties and 
we describe the tools and methods for their analysis. In sec- 
tion V we briefly describe the Fermi-LAT instrument and 
the characteristics relevant for our analysis. In section VI we 
discuss the results of our data simulation pipeline. Finally, 
in section VII we evaluate Fermi's sensitivity to distinguish 
a scenario in which the diffuse gamma-ray emission has a 
DM contribution (either extragalactic or galactic in origin) 
with respect to the case in which DM is absent. Conclusions 
are presented in section VIII. 



2 GAMMA-RAY SKY COMPONENTS 

2.1 Isotropic Gamma- Ray Background 

In the following, for clarity, we will refer to the isotropic ex- 
perimentally detected emission as the Isotropic Gamma-Ray 
Background (IGRB) given that in principle its origin may 
be not fully extra-galactic. Indeed, we will also consider the 
possibility of a nearly isotropic emission from DM annihila- 
tion in the Milky Way substructures. We will use the term 
Extra-Galactic Gamma Background (EGB) when referring 
to predictions from models of extra-galactic emission either 
from DM annihilation or from astrophysical processes. 

A first detection of the IGRB in the GeV range was re- 
ported by the EGRET experiment ([Sreekumar et al. 1998)1 
(see also the analysis by Strong, Moskalenko and Reimer 
( [Strong et al. 2004a[ ) ). The recent measurement from The 
Fermi collaboration ([Abdo et al. 2009b)) finds that the 
IGRB energy spectrum has an almost perfect power law 
with slope s ~ —2.4 which is approximately what is expected 
from the contribution of unresolved blazars thus suggesting 
an extragalactic origin of the IGRB. 

A (possibly incomplete) list of contributions to 
the IGRB is given by the emission in astrophysi- 
cal processes from Blazars, AGNs, Starburst Galax- 
ies, Star Forming Galaxies, Galaxy clusters. Clusters 
Shocks, Gamma-Ray Bursts ([Sleeker and Salamon 19961 
IPavlidou and Fields 20021 IGabici and Blasi 20021 

ITotani 1998p . A contribution is also expected from ultra 
high energy photons and electrons cascading to GeV-TeV 
energies through interactions with the CMB or the Infrared- 
Optical Background ([Kalashev et al. 2007)) . Each of these 
source classes is expected to produce anisotropies although 
all them should more or less trace the filamentary pattern 
of the cosmological Large Scale Structures. Anisotropies 
will be discussed in more detail in the next two sections. 



2.2 DM Annihilation 

Apart from the astrophysical contributions, an exotic com- 
ponent arising from WIMP annihilation could also be ex- 
pected. In the case of uniformly distributed DM, the an- 
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Figure 1. Energy spectrum of the various components contributing to the gamma sky above 30° (left) and above 60° (right) galactic 
latitude. The isotropic component is the determination from Abdo et al. (2009b). Also showrn is an example of extragalactic WIMP 
emission normalized with the model of UUio ct al. (2002) (with NFW halo profile and no halo substructures). The spectrum of the 
resolved point sources is calculated from the Fermi 11-month source catalogue (Ballet 2009). 



nihilation flux for a thermal relic cross section ( (aAv) ~ 
3 X 10~'^®cm"^s^^ ) is many orders magnitude below the 
IGRB. However, given that the DM distribution on cos- 
mological scales could be very clumpy, the enhancement 
factor (denoted as A'^{z)) due to the quadratic depen- 
dence of the annihilation signal on the DM density has 
to be taken into account. This enhancement is typically 
several orders of magnitude, thus boosting the DM signal 
to a level comparable to the IGRB ([Bergstrom et al. 2001[ 



lUUio et al. 2002l[Tapor and Silk 2002[ ). The exact boost de 
pends on the details of the modeling of the DM clustering, 
the DM halo profile and the structure formation history. 
We fix the absolute normalization of the DM spectra fol- 
lowing the DM haloes clustering model described in UUio et 
al.(2002). Two versions of the model are employed: a con- 
servative version assuming a Navarro-Frenk- White (NFW) 
profile (|Navarro et al. 1995|) for the haloes and no substruc- 
tures, which gives a low DM normalization, and a more op- 
timistic version where the haloes still follow a NFW profile 
but the DM signal is enhanced by the presence of substruc- 
tures in the DM haloes. It should be noted that the clus- 
tering properties of the DM haloes are still very uncertain, 
especially at small scales (approximately below the typical 
galactic halo scale), and even larger normalizations are pos- 
sible. For example recent results from the Millenium-II sim- 
ulation (|Zavala et al. 2009|l indicate that, depending on the 
extrapolation employed for the low mass haloes (non re- 
solved in the simulation), an order of magnitude enhance- 
ment with respect to the above "optimistic" case can be 
achieved. Overall, Zavala et al. (2009) place the possible val- 
ues of the boost A'^(z=l) in the range 10* to 10^ (see also 
(|Abdo et al. 2010a|) . in particular Fig.l). 

2.3 Galactic Foregrounds 

The main contribution to the Galactic foregrounds is given 
by the decay of pious produced in the interaction of Cosmic 



Ray (CR) nucleons, by the inverse Compton (IC) scattering 
of CR electrons on the interstellar radiation field (ISRF) and 
by bremsstrahlung of CR electrons on the interstellar gas. A 
further contribution from pulsars has also been considered 
recently ( [Faucher-Giguere fc Loeb 2009^ . According to this 
analysis, the contribution from Galactic millisecond pulsars 
can be quite isotropic (i.e. can contribute to high Galac- 
tic latitudes) and constitute a large fraction of the IGRB. 
However, lacking more precise studies of the subject, we will 
neglect this contribution for the present analysis. 

The normalization of the Galactic foregrounds is 
clearly important to identify the regions where the 
IGRB is dominant, i.e. the regions most promis- 
ing for the anisotropy analysis. We employ the so 
called "conventional model" ( [Strong and Moskalenko 1998[ 
[Moskalenko and Strong 1997[]Strong et al. 2004b[ ) which is 
derived under the assumption that the CR electron and nu- 
cleon spectra in the galaxy can be normalized by the lo- 
cal (solar system) measurements. This model represents a 
nice fit of the Fermi data, at least at Galactic latitudes 
|6| > 10°, where the emission is mostly local in origin (com- 
ing within a few kpc's of the solar position) (|Porter 2009al 
I Abdo et al. 2009cl [Porter 2009b|) . We wiU, indeed, focus on 
this region of the sky in the following. The Galactic center re- 
gion and the Galactic plane are anyway strongly dominated 
by the Galactic emission and they are thus not suitable for 
the anisotropy analysis. A description of the diffuse emission 
in terms of the conventional model is thus accurate enough 
for the present analysis. 



2.4 Resolved Point Sources 

The number of resolved point sources has grown consider- 
ably in the Fermi era relative to the roughly 270 in the 3rd 
EGRET catalogue (IHartman et al.(1999)l). The first year 



Fermi catalogue includes 1451 sources with a significance 
larger than about 4 sigma (|Abdo et al. 2009d)) . We include 
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explicitly these sources in our analysis. As will be clear in 
the following, it is important to include and, in addition, 
to model this component correctly since it is a considerable 
source of anisotropy which needs to be reduced to gain sen- 
sitivity to the anisotropies of the diffuse component. 

While we use the first year catalogue, in the following we 
produce a forecast over 5 years of data taking, thus should 
in principle be using an associated 5 years catalogue. The 
effect of using a five year catalogue, however, is expected to 
be small: the masked area will increase, but on the other 
hand the Poisson noise will be reduced. Both effects are 
small since increasing the data taking time for the catalogue 
will add sources on the faint end, requiring small masking 
area and producing a small decrease in Poisson noise (since 
the brightest sources are already removed). On the other 
hand, modeling the 5 year catalogue requires an assumption 
on the nature of the unresolved sources, introducing a model 
uncertainty. We therefore decided to use the 1 year catalogue 
only. 

2.5 Unresolved Point Sources 

Part of the IGRB is made of unresolved point sources. It is 
still a matter of debate if unresolved point sources can make 
up all the IGRB or only part of it. This issue is likely to be- 
come clearer when population studies of the Fermi resolved 
point sources become available. Present studies indicate that 
unresolved blazars can only account for at most 30 % of the 
IGRB measured by the Fermi- LAT (|Abdo et al. 2010b|l . 

Unresolved point sources can also be a further source of 
anisotropy in the IGRB itself, producing a Poisson noise-like 
contribution to the anisotropy spectrum. It is worthwhile to 
point out that this Poisson noise is different from the Poisson 
shot noise in the angular power spectrum (see below and 
appendix A), which is related to the finite number of counts. 
The contribution of unresolved point sources depends both 
on the ability of the Fermi-LAT to detect point sources and 
remove them from the diffuse emission and on the intrinsic 
number density of the sources themselves. Both these terms, 
in turn, are quite model dependent but they are expected to 
be reasonably constrained from population studies of source 
classes with resolved members. 

For the moment, for the sake of studying the impor- 
tance of the effect, we consider, beside the model described 
in the next section, a model fully made of unresolved point 
sources. For the Poisson noise we assume a value of Cp = 
5.0 X 10~^ which is in the middle between the estimate of 
10"* (|Ando et al. 2007|) (see their Fig.4) and the prelim- 
inary measurement of 10^'' from the Fermi collaboration 
( jSiegal-Gaskins 2010[ ). 

It is also worth noting that the Poisson noise of unre- 
solved sources is a time varying term slowly decreasing in 
time. The brighter unresolved sources, in fact, get finally 
resolved, as Fermi continues to collect more statistics, and 
once they are removed, and/or masked, the "new" remain- 
ing unresolved signal has a slightly lower level of Poisson 
noise. A detailed study of this effect requires extrapolation 
from population studies which we leave for future work. 

In Figure [l] we show the energy spectrum of the various 
Galactic and extra-galactic components for the model we 
employ. In particular, we show the average spectrum in the 
regions |&| > 30° and |b| > 60°, where the IGRB has the 



highest relative contribution. We show the power law fit to 
of the IGRB from (|Abdo et al. 2009b|) . Also shown is an 
example of extragalactic WIMP emission normalized to the 
model of Ulho et al. (2002) (with NFW halo profile and no 
halo substructures). The contribution from resolved point 
sources detected by Fermi (jAbdo et al. 2009a|) . (|Ballet 2009|l 
is also present. It can be seen that, combined, these sources 
constitute a significant fraction of the diffuse background. 
For simplicity, the spectrum of each single source is modeled 
as a simple power law and thus the behavior at high energies 
( > 10 GeV) may have some biases with respect to the real 
spectrum. 



3 ANISOTROPIES: GENERAL 
CONSIDERATIONS 

We will follow two approaches to quantify the anisotropies 
of all the gamma-ray sky components. For the components 
for which maps as a function of energy are available (i.e. for 
the Galactic foregrounds and for the resolved point sources) 
we calculate the angular power spectra directly from the 
maps using HEALpix (|Gorski et al. 2005|l (see appendix A 
for more details on HEALpix and the definition of power 
spectra). For the remaining components (DM and astrophys- 
ical EGB) we model the power spectra from theory. From 
the theoretical power spectra we then also produce simu- 
lated sky-maps (under the assumption that the fluctuations 
are gaussian). In this way we have effective sky maps for all 
the components and we can combine them, convolve them 
with the Fermi IRFs, simulate finite statistics of events, ap- 
ply masks, etc. The final maps will thus represent a sim- 
ulated gamma-ray sky as seen through the Fermi detector 
and where all the components are entangled. Power spec- 
tra are then recovered at the end from the final maps using 
HEALpix. In this way the sensitivity to DM can be studied 
as described in section [7] by comparing the power spectra of 
the maps with and without the DM component. 

Fig[2] shows separately the all-sky angular spectra of 
the various components: Galactic foregrounds, astrophysical 
EGB, WIMPs EGB and point sources. The spectra are cal- 
culated from simulated maps at _E > 1 GeV after convolving 
with the Fermi IRFs (see the next sections). For illustrative 
purposes, the flux of the DM component is arbitrarily nor- 
malized to the level of the IGRB detected by Fermi and an 
observation time of 5 years is assumed. It can be seen that 
the Galactic foreground anisotropies have a decreasing spec- 
trum in I and dominate over the extragalactic anisotropies 
(either astrophysical or WIMP) up to i ~ 200 - 300. Luck- 
ily, however, this can be remedied using an appropriate mask 
excluding the Galactic plane region where most of the fore- 
grounds are concentrated. In this way the foregrounds can 
be reduced such as to become subdominant at a multipole 
of i ~ 30 as shown in the lowest panel of figH] The spec- 
trum of point sources highly dominates both the foregrounds 
and the extragalactic signals. This indicates that an efficient 
point source masking is crucial to extract the interesting ex- 
tragalactic signal. These steps will be explored in the next 
sections. The figure also illustrates the effect of the Poisson 
noise (first panel) coming from finite statistics present in the 
calculated raw power spectra. Customary power spectra are 
plotted with the noise removed (second panel). 
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Figure 2. All sky angular spectra for E > 1 GeV of the employed 
models for point sources (rod), Galactic foregrounds (green), EGB 
(magenta) and EGB DM (blue) after convolution with the Fermi- 
LAT Point Spread Function. For illustration, the flux of the DM 
component is arbitrarily normalized to the level of the IGRB de- 
tected by Fermi- LAT, and 5 years of Fermi-LAT observations have 
been assumed. The power spectrum of the exposure map is also 
shown (black line). In the top panel the spectra before shot noise 
removal are shown (shot noise is represented by the dashed lines) 
. All the spectra are calculated with HEALpix. The last panel 
shows the angular power spectra after the application of a suit- 
able mask to cover the low latitude Galactic foregrounds and the 
point sources. Since the mask is effective in suppressing the point 
sources signal the related power spectrum has been removed. 



The following section describes in detail how the power 
spectra of DM and astrophysical EGB are modeled and how 
the related maps are simulated. The section can be skipped 
by the reader not interested in these details. 



4 DARK MATTER AND ASTROPHYSICAL 
ANISOTROPIES 

4.1 Modeling the EGB 

Apart the Poisson term coming from the unresolved point 
sources, the remaining source of anisotropics of the IGRB is 
given by the anisotropic spatial distribution of the sources 
themselves. To derive the anisotropy we will assume, as a 
reasonable first approximation, that the gamma ray sources 
are distributed as the matter density of the universe p{x), 
i.e. following the cosmological Large Scale Structures (LSS). 
In principle ps, the density distribution of astrophysical 
sources, should be used instead of p: ps in general exhibits 
a scale and time dependent bias with respect to the matter 
density. However, specific classes of astrophysical gamma- 
ray sources have different biases. For example, blazars are 
well known to concentrate at the center of clusters of galax- 
ies, thus presenting an over-bias with respect to galaxies at 
high densities. On the other hand, galaxies and clusters of 
galaxies reasonably trace the matter density, at least in the 
recent cosmic epoch. The assumption p^ = p is thus general 
enough to approximately describe emission from astrophys- 
ical sources. 

Given these assumptions the extragalactic cosmic 
gamma-ray signal can be written as (|Ullio et al. 2002] 
[Bergstrom et al. 2001) ICuoco et al. 2006|l 



I-y(E.y, h) CX 



p{z,h,r{z))g[E,{l + z)]e-^('=^'^'> 
H{z){l + zy-' 



(1) 



where g{E) — dN-,/dE is the photon spectrum of the 
sources, E^ is the energy we observe today, p{z,n,r) is the 
matter density in the direction n at a comoving distance 
r, and the redshift z is used as time variable. The Hub- 
ble expansion rate is related to its present z — value Ho 
through the matter an d cosmological con stant energy den- 
sities as H{z) = //o\/f2i\/(l -|- z)^ + Q,A, and the reduced 
Hubble expansion rate h{z) is given by H{z) — 100 h{z) 
km/s/Mpc. We will in the following use Qm — 0.3, Q,a ~ 0.7 
and Ho = 70 km/s/Mpc. The quantity t{E^,z) is the op- 
tical depth of photons to absorptions via pair production 
(PP) on the Extra-galactic Background Light (EBL). We use 
the parametrization of T{E-y,z) from (|Stecker et al. 2006|l 
for < 2 < 5, where the evolution of the EBL is included 
in the calculation. The EBL is expected to be negligible at 
redshifts larger than 2 ~ 5 corresponding to the peak of 
star formation. Thus, gamma photons produced at earlier 
times experience an undisturbed propagation until z ~ 5, 
while only in the recent epoch they start to lose energy 
due to scattering on the EBL. Correspondingly, we assume 
t[E^,z) = T{E.y,5) for z > 5 (see also formula (A. 6) in 
(|Cuoco et al. 2006|) ). 

In the case of cosmological DM annihilation, the re- 
sulting spatial distribution of the gamma signal follows the 
square of the matter distribution p'^{x) through 
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IxiE^,n) 



<axV> 



pl{z,h,r{z))g[E-,{l 



-T(Ey,z) 



H(z){i + zy 



,(2)Cx 



dr 



Kir) P, 



{k=L,zir)) 



(8) 



which, taking the spatial average, reduces to the expression 



- 2 roo ,2 



IxiE-,) = 



<iTxV> pa 
87rm| jQ 



H{z) 



where A^(2;) is the dumpiness enhancement factor and po 
is the average DM density at z = 0. 

We can write the astrophysical EGB and DM in a com- 
pact form: 

/•oo 

I^{Ei,E2,n) oc / zW^{Ei,E2,z)p{z,h), (4) 

Jo 



(5) 



Ix{Ei,E2,h) oc / zWx{Ei,E2,z)p\z,n) 



where W.y{Ei,E2,z) and Wx{Ei,E2,z) are the astrophysi- 
cal and DM window functions, which contain the informa- 
tion about gamma-ray propagation, injection spectra and 
cosmological effects 



W{EuE2,z) 



l' g[E{l + z)]{l + zf-\ _^^^^^^ 

:• H{^) 



(6) 



where a — 1,2 applies in the astrophysical and DM cases, 
respectively and Ei and E2 are the boundaries of the energy 
bands considered. We are using the notation p{z,n,r{z)) = 
[l+z)'^ X p{z, n) to underline that the window function is only 
dependent on the two variables, direction and redshift, and 
to make the (l+z)"^ behavior of the matter density explicit. 
In particular, choosing Ei — Ecut and E2 = -|-oo, when 
properly normalized W{Ecut,z) represents the probability 
of receiving a photon of E-, > _Ecut emitted at a redshift z. 
This can be used to define an effective horizon, z-h, beyond 
which the probability of receiving a photon is negligible. In 
practice, above about 10 GeV, pair production attenuation 
sets the scale of the horizon which is about zn ~ 0.5 and 
z-H ~ 2 for the a = 1,2 cases respectively. In this case, 
DM annihilation produces a much larger anisotropy than the 
astrophysical sources. In contrast, below 10 GeV the horizon 
is set by cosmological effects and z-h ~ 1 and z-h ~ 10 for 
the two cases. In the low energy regime we thus see that the 
DM anisotropics are averaged over a much larger horizon 
and the final signal can be smoother than the astrophysical 
case. Further details are discussed in (| Cuoco et al. 2007)l . 

The injection energy spectrum for DM annihila- 
tion g{E^) = dN^/dE, given the annihilation chan- 



SUSY (|Gondolo et al. 2004|) . which, in turn, uses a tabu- 
lation derived from simulation of the particle processes per- 
formed with Pythia ( |Sjostrand et al. 2007[ |. We will consider 
in the following the bb, t'^t~ and p^ pT annihilation chan- 
nels. 



4.2 Angular Power Spectra 

Given the above expressions we can write the angular 
power spectra of the various dimensionless fluctuation fields 
5II<I> 



C\ = j%W',{r)P,(k=i,z{r)) 



(7) 



where we have used the Limber approximation 
N (|Limber 1953[l . which is accurate for all but the very 
, (tawest multipoles, while Pp{k,z{r)) and Pp2{k,z{r)) are 
3D power spectra of the matter field and of its square. We 
are also interested in the spectra of cross correlation among 
two different energy bands, which can be written as 



C12 



% Wi{r)W2{r) 



P k 



,z(r) 



(9) 



where i — p, p^ for the astrophysical and DM case. 

We take the 3D power spectra Pp and Pp2 from an 
N-body simulation of Large Scale Structures formation 
(|Cuoco et al. 2007)1 . We note that Pp2 has a dependence 
on the clustering properties of DM below galactic scales 
(see (|Ando fc Komatsu 2006 |l for an analytical derivation 
of Pp2) which are not resolved in the N-body simulation. 
The anisotropy properties in the case of cosmological DM 
annihilation with substructures are therefore not fully con- 
sistent with its energy spectrum normalization, which does 
take into account the effects of substructures. This has to be 
taken into account when interpreting the sensitivities in the 
cosmological substructures model shown in the following. 

As described in the previous section we will consider 
in the following two different anisotropy models for the as- 
trophysical EGB in order to bracket our uncertainties and 
estimate the effect on the DM sensitivity. The first one is a 
model dominated by the correlation term of the power spec- 
trum (EqO, while the second one is a model dominated by 
the Poisson term Cp from unresolved point sources. In this 
case the angular power spectrum is simply: 



C-y = Cp 



(10) 



with Cp = 5 X 10"^ 

Finally we will also consider the possibility that 
part of the IGRB is produced by DM annihilation in the 
sub-structures of our galaxy as found by recent N-body 
simulations (jDiemand et al. 20071 [Springel et al. 2008a| 
[Springel et al. 2008bp . In this case, compared to the cosmo- 
logical annihilation scenario, the situation is simpler since, 
in the absence of redshift effects, the pattern of anisotropics 
is independent of energy (and thus of the particle physics 
model) . This can be seen more clearly in Eq[S] where W 
becomes a constant independent of energy and redshift. The 
anisotropics thus only depend on the clustering properties 
of DM in the Milky Way (MW). We will take as a reference 
model the C^gai from ( [Siegal-Gaskins 2008[ ) where the 
anisotropics have been calculated from a simulated MW 
substructure map. In particular, we will use the spectrum 
shown in their Fig. 5 for a minimum clump mass of IOM0. 
Compared with the cosmological DM anisotropics this 
model gives a normalization of the angular power spectrum 
a factor ~ 100 larger (or more) depending on the energy. 
More precisely, the shapes of the power spectra are sim- 
ilar for the two cases, while the absolute normalization 
for the Galactic case is l{l + 1)Ci/2tt « 10 at Z = 100 
independent of energy, whereas for the extragalactic case 
l{l + l)Ci/2n « 0.02 at 10 GeV. Fig|3] shows the energy 
dependence of l{l + l)Ci/2n at I = 100 for the various cases 
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Figure 3. Energy dependence of the anisotropies for our models 
of cosniological EGB DM, astrophysical EGB with and without 
a dominant Poisson term and DM in Galactic substructures. The 
anisotropy curve for the DM refers to a m^=l TeV WIMP anni- 
hilating into bb. 



while FigHshows l{l + l)Ci/2n for the DM EGB case for 4 
difTerent energies. 



i,j = 1,Ne for the DM or astrophysical EGB. We then 
generate the maps in harmonic space («;„, (see definitions 
in the appendix)) up to Z ~ 500 sampling from the 60x60 
matrix of the multi-gaussian distribution given by the C,'-' 
for each I. In practice the energy behavior is generally highly 
correlated so that only the first 4-5 principal components are 
required for the simulation of the 60 maps. The synf ast tool 
from HEALpix is then employed to combine the harmonic 
components into the real space maps /b. (f2). These energy 
dependent anisotropy maps are then weighted and normal- 
ized at each energy according to the energy spectrum of the 
component itself (astrophysical or DM EGB), giving the fi- 
nal template maps. These can now be added to the model 
of Galactic foregrounds and to the model of point sources 
to give a simulation of the gamma-ray sky where all the 
components are entangled. 

For the case of Galactic DM the procedure is simpler 
given that the anisotropy pattern is independent of the en- 
ergy and only one map needs to be simulated. The same 
holds for the astrophysical EGB model with pure Poisson 
noise. 

Convolution with the instrumental effects is described 
in the next section. 



4.3 Model Sky Maps 

To produce a realistic simulation of data from Fermi 
we have to simulate the expected sky pattern of every 
component. In the case of Galactic foregrounds we em- 
ploy the GALPROP package ( [Strong and Moskalenko 1998[ 
[Moskalenko and Strong 1997[ ) which produces as output the 
sky maps of the foregrounds as a function of energy. As 
already described above the conventional diffuse model is 
used. It is worth noticing that energy dependent CR diffu- 
sion effects imply an energy dependent sky pattern as well. 
In general, at low energies the diffusion length of CRs is 
large and the resulting maps are thus smoother than the 
higher energy counterparts where the diffusion length de- 
creases. Thus the relative importance of the Galactic fore- 
grounds anisotropies at small scales {I > 100) increases with 
energy. Of course, projection effects along the line of sight 
can sometimes complicate the picture and give small scale 
anisotropies also at low energies although this in practice 
happens only along the Galactic plane where there is a con- 
tribution from structures many kpc away from the solar po- 
sition. 

For the remaining components, i.e. the astrophysical 
EGB and the DM EGB, we create template maps as gaus- 
sian realizations of the theoretical angular power spectra C' 
derived above. The situation is however complicated by the 
fact that the angular spectra are energy dependent and, fur- 
ther, that there are non-zero cross-correlation power spec- 
tra between different energy bands. Energy dependent maps 
with the correct cross-correlation properties among differ- 
ent energies need thus to be realized. In practice, we choose 
Ne = 60 logarithmically spaced energy bands between 10 
MeV and 10 TeV (i.e 10 bands per energy decade). We 
then calculate the angular auto-correlation power spectra 
for each band and the cross-correlation power spectra among 
all the possible pairs of energy bands according to equations 
I7I10I The result is the model power spectrum matrix C^'' 



5 THE FERMI-LAT AND SIMULATED 
DATASETS 

In the following we describe the various instrumental effects 
introduced in the Gamma-ray sky when observed through 
the eyes of the Fermi-LAT telescope. 



5.1 Charged particles contamination 

Charged particles, mainly protons, electrons and positrons, 
as well as a smaller number of neutrons and Earth albedo 
photons, present a major instrumental background to poten- 
tial DM signals, in particular when considering the isotropic 
gamma-ray signal. These background particles greatly dom- 
inate the flux of cosmic photons incident on the LAT, and a 
multivariate technique, employing information from all LAT 
detector systems, is used to reduce the remaining charged 
particle background to only a fraction of the expected 
flux of extragalactic diffuse emission (|Atwood et al. 2009p . 
The analysis of the extragalactic background from Fermi 
(|Abdo et al. 2009b)) however indicates that, with the stan- 
dard data selection cuts, the residual contamination is com- 
parable to the photon background itself above a few GeV 
and dominates at higher energies. Generally, this back- 
ground is expected to be basically isotropic so it can be easily 
taken into account in the anisotropy analysis. Alternatively, 
stringent cuts can be applied on the photon events decreas- 
ing the contamination but also somewhat the exposure and 
thus the signal (|Abdo et al. 2009b)) . For our simulation we 
will use the effective area for the "diffuse" event class as de- 
fined in (|Atwood et al. 2009|l . Using more stringent selection 
cuts to reduce the effects of CR contamination reduces the 
effective area by a factor of ~25% (|Abdo et al. 2009bp . The 
sensitivity curves should change accordingly, approximately 
by the small factor Vl-25 ~ 1.1. This estimate is reasonably 
accurate for a sensitivity study as the present one. We do 
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not attempt a detailed simulation of the CR background it- 
self, postponing a more careful investigation of the issue to 
the data analysis work. 



5.2 Point Spread Function and Exposure 

The effects of the Fermi-LAT detector can be summarized 
in terms of the IRFs which include the point spread func- 
tion (PSF) and the energy dependent exposure maps. For 
this analysis we use the PSF and exposure as available from 
the Fermi Science tools p6v3 Q We neglect the effect of en- 
ergy dispersion (which is about 10 %) and assume that the 
photon energy is perfectly reconstructed. In this respect, for 
this analysis we will divide the simulated datasets in wide 
energy bands so that this instrumental effect is negligible 
for our purposes. To obtain simulated maps we use the en- 
ergy dependent model input maps and then convolve them 
with the Fermi PSF and simulated exposure of 5 years. More 
specifically, the exposure was generated from 1 yr of pointing 
history simulated with gtobssim and then rescaled to 5 yrs. 
The convolution with the PSF and the exposure was then 
performed using the GaDGET (|Ackermann et al. 2009[) tool 
(see also Abdo et al. (2009b)). The resulting maps have 
a normalization in units of counts per steradian which is 
further re-normalized to counts per pixel (~ 10~^sr for a 
Nside=256 HEALpix pixelization). Finally Poisson random 
noise is added to simulate the effect of finite statistics. 

The Fermi PSF is energy dependent and improves at 
high energies. It corresponds to an angular resolution of 
about 3° at 100 MeV, 0.6° at 1 GeV and 0.1° above 10 GeV 
(see (|Atwood et al. 2009|) for more details). The PSF thus 
heavily affects the anisotropies especially below 1 GeV, sup- 
pressing the medium-high multipoles in the angular power 
spectrum. The sensitivity to these multipoles is not com- 
pletely lost however. In principle, from the knowledge of the 
PSF it is possible to recover the true power spectrum even 
at higher multipoles although the error grows with larger 
attenuation (it grows exponentially with multipole I: the 
suppression factor is given by Wi — exp{—l^a^/2), for a 
gaussian PSF with variance at )■ 

The effect of the PSF can be seen in Fig|4] where our 
model power spectra for the DM EGB component at vari- 
ous energies are compared with the corresponding PSF con- 
volved spectra. The latter are obtained by convolving the 
model DM EGB sky maps (simulated from the angular spec- 
tra with the procedure described in the previous section) 
with the PSF and extracting the power spectra from the 
resulting maps. As expected, the PSF heavily suppresses 
anisotropies at I > 100 for E < IGeV. Notice also that there 
is a certain suppression at I > 200 even for E — 300 GeV, 
which is slightly more than the factor Wi = exp{—l^a1/2) 
expected at this energy for a gaussian PSF with width 0.1°. 
This extra effect is likely related to the non-gaussian form 
of the PSF which even at very high energies (> 100 GeV) 
exhibits a tail towards large angular spreads. The figure also 
shows the difference in PSF for events converting in the front 
or the back of the detector where the reconstruction perfor- 
mances are different. It can be seen that the PSF is better 
for events converting in the front, giving a sensitivity to 



See http://fermi.gsfc.nasa.gov/ssc/ 




Figure 4. Power spectra at various energies of the model of DM 
EGB emission (for a m^=l TeV WIMP annihilating into bb) and 
their corresponding PSF convolved counterparts (from top to bot- 
tom 3.t E = 300 GeV, B = 30 GeV, E = 2, GeV, E = 300 MeV). 
The PSF suppresses anisotropies at I > 100 for E < IGeV. The 
rise of the power spectrum at high multipoles in the E = 300 MeV 
case is a numerical HEALpix artifact. The lower panel shows the 
power spectra resulting from the different PSFs for events con- 
verting in the front or the back of the detector and for the average 
PSF. The dipole has not been included in the simulated maps, 
hence the sharp rise at I = 2. 



slightly higher multipoles especially at low energies. Above 
a few GeV, instead, the difference is less important. We will 
use the average PSF for our analysis. 

The exposure maps are very uniform over an averaging 
period of 1 year (see again (|Atwood et al. 2009p ) so that the 
convolution with the exposure only adds power to the very 
lowest multipoles I < 10. This can seen in fig[2]where we plot 
the angular power spectrum of the exposure which, indeed, 
falls rapidly at / = 10. The exposure pattern is also energy 
dependent, although only weakly. The integrated exposure, 
however, rises steeply starting from about 100 MeV and then 
fiattens above about 1 GeV. 



6 RESULTS OF THE SIMULATION 

6.1 Simulated Sky-Maps 

The maps produced as described in the previous sections are 
summed up to build a simulation of the observed gamma- 
ray sky. A reference no-DM model is derived adding Galactic 
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Figure 5. (Left) 5 year averaged counts map (in HEALpix nside=256 format) of our reference no-DM model, which contains astrophysical 
EGB, Galactic foreground and resolved point sources. An energy dependent mask is applied to suppress Galactic foregrounds and point 
sources (see the text). The simulation has been divided into five energy bands (in the energy ranges indicated in the titles) and only the 
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E = [0.4,0.8] GeV 
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Figure 6. Power spectra matrix for the 5 average counts maps (times the maslcs) of the previous figure. Before computing the spectra, 
the maps have been (re)normalized so that to represent fluctuation maps. Auto-correlations are in the diagonal and cross-correlations 
are off-diagonals. The dashed lines in the auto-correlation plots indicate the level of shot noise oc l/A^counts- The noise in different energy 
bands, in contrast, is uncorrelated so that no shot noise is present in the cross-correlations. The blue lines represent the power spectrum 
matrix of the masks alone. It can be seen, that the masks give the "bulk" level of anisotropics. The anisotropics of the masked no-DM 
reference model then add over the mask level (black curves). At higher energies the angular resolution of Fermi-LAT increases and the 
anisotropy excess with respect to the masks is better visible (in the lower right part of the plot) at medium-high multipoles. To be at a 
significant level, the difference has to exceed the errors bars. Statistical errors are represented by the gray shaded area for a confidence 
level of two sigma. Finally, the red curves represent the angular power spectra when a DM EGB component is added on top of the 
previous contributions. For illustration m-^ ft! 200 GeV is chosen and the DM flux is normalized to the same level as the astrophysical 
EGB at 10 GeV. The bump visible above the statistical error level in the 10 GeV-100 GcV autocorrelation spectrum (and the absence 
of it in the other energy bins) would make this model detectable. 
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foregrounds, astrophysical EGB and point sources. The as- 
trophysical EGB is normalized to the level of the observed 
IGRB. The normalization of the Galactic foregrounds is the 
one given as output by GALPROP. Point sources are nor- 
malized to their observed flux. The overall statistics is nor- 
malized to 5 years of data taking by Fermi. The maps we 
obtain are average counts maps, but a realization of effec- 
tive counts maps can be easily obtained drawing poisson 
counts in each pixel from the average counts maps. Although 
we use 60 energy bands to produce the simulated maps, it 
is convenient to rebin them into fewer maps, in order to 
have a reasonable amount of statistics in each band. We 
thus divide the dataset in 5 energy bands from 400 MeV 
to 800 GeV. Finally, a crucial step is the construction of 
the masks which will suppress the Galactic foregrounds and 
the resolved point sources. For each energy band the mask 
are defined as the region where the foregrounds do not ex- 
ceed twice the EGB (200%), and the point sources emission 
does not exceed more than 20% of the EGB itself. Simu- 
lated sky maps (outside the masked regions) for both the 
average counts and one particular realization of the counts 
are shown in Fig[5] Indeed, outside the masks the remaining 
signal is quite isotropic, indicating that the IGRB makes 
most of the signal, but still a certain residual foreground 
contamination is apparent, especially near the edges of the 
mask. Clearly this is expected given the loose cut employed 
to build the masks (only 200% of the EGB). A tighter cut 
on the foregrounds can be applied but then the available 
area for the analysis would shrink to close to zero since the 
Galactic foregrounds never go below 20-30% of the EGB 
even at the Galactic poles. For the present forecast anal- 
ysis a substantial foreground contamination is still accept- 
able since the foreground model is known in advance and 
what is required for a significant detection is that the EGB 
anisotropy exceeds the anisotropy of the residual Galactic 
foregrounds. In general, for a real data analysis some kind 
of foreground cleaning must be performed to obtain useful 
results. We postpone a detailed treatment of this issue to 
the analysis of the real data. 



On the other hand, the cut of 20% of the EGB for point 
sources is rather effective and removes most of the gamma- 
ray events from resolved point sources, except perhaps few 
events around the most powerful sources where, statistically, 
some events very far from the source position are expected 
due to the tails of the PSF. From Fig[5]it can also be seen 
that the fraction of sky where the isotropic component domi- 
nates increases with energy. This is expected considering the 
fact that the IGRB has a spectrum harder than the Galac- 
tic emission (see also Fig[T]). Also it is clear that the area 
masked for the point sources substantially decreases at high 
energy due to the improvement of the PSF. This is true for 
all but the highest energy band, where the fraction of sky 
covered by the point sources increases again. This can be 
understood looking at the energy spectrum of the resolved 
sources in Fig[l] This becomes very hard for -B > 10 GeV 
so that, even if the PSF at this energies is very narrow, its 
tail still contains a flux which is a significant fraction of the 
IGRB. 



6.2 Angular Power Spectra 

In Fig[6] we plot the various auto-correlation and cross- 
correlation spectra in a "triangle plot", where the diagonal 
represents the auto-correlations and the the off-diagonals 
the cross-correlations. The blue lines represent the spectra of 
the masks only, which thus provide the zero-level anisotropy. 
The black line represents the average spectrum of the ref- 
erence no-DM model convolved with masks, i.e. the maps 
in the left column in Fig[S] As explained above, outside the 
mask the reference model is to a large extent given by the as- 
trophysical EGB, although a substantial residual contamina- 
tion from the Galactic foregrounds is still present. The resid- 
ual point source contamination, on the other hand, is neg- 
ligible. The grey shaded area is the two sigma error region 
in the power spectrum accounting for the effect of the finite 
number of events. A logarithmic binning of the multipoles 
in 18 bands is used. The level of shot noise is shown with a 
dashed line. It can be seen that the shot noise increases in 
the high energy bands where the statistics decrease, as ex- 
pected. Finally, the red curves represent the angular power 
spectra if a EGB DM component is added on top of the pre- 
vious contributions. For illustration a WIMP with m^ ~ 200 
GeV is chosen and the DM flux is normalized to the same 
level of the EGB at 10 GeV. It can be seen, indeed, that in 
the 4th energy band (approximately 10-100 GeV), where the 
DM emission peaks, the anisotropy has a noticeable bump 
compared to the EGB anisotropy and provides a clear sig- 
nature for this model. The DM anisotropy signal disappears 
in the 5th band (100-800 GeV) where the DM contribution 
drops. This "energy modulation in the anisotropy spectrum" 
is also discussed in ( [Hensley et al. 2009[ )). 

It has to be noted that the spectra look quite different 
from the all-sky spectra of the astrophysical and DM EGB 
shown in Fig[2] both in the shape and normalization. This 
is of course due to the effects of the residual foregrounds 
and especially of the mask which distorts the original spec- 
trum. The mask effects can in principle be removed with 
specifically designed algorithms like the MASTER algorithm 
(|Hivon et al. 200"2)) . This comes, however, at the price of fur- 
ther errors in the final spectrum which need to be addressed 
with Montecarlo simulations (a package implementing the 
full procedure is described for example in (|Lewis 2008|) ). 
For the present purpose, it suffices to use the masked sky 
power spectrum without attempting to deconvolve it from 
the mask itself. This is fine as long as all the power spectra 
are calculated on the same sky region in order to allow a 
consistent treatment. 



7 DM SENSITIVITY 

It is in principle straightforward to calculate the sensitivity 
to DM requiring that the DM anisotropy does not exceed 
the anisotropies of the reference no-DM model. Let us de- 
fine C^Qgj, as the binned multipoles power auto- and cross- 
correlation spectra of the reference model Q The indices i,j 
run from 1 to 5 and indicate the energy bands and the index 



■^ We label it for simplicity EGB although, besides the true as- 
trophysical EGB, there are some residual contaminations from 
Galactic foregrounds and point sources (see previous section). 
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Figure 7. Exclusion plots in the my^-{crp^v) plane for various 
annihilation channels: bb, t'^t~ , fj,'^ ^~ . Each panel shows the 
sensitivity curves at the 2(t, 3o" and Scr levels. Two different nor- 
malizations for the flux of the cosmological DM signal have been 
employed using an updated version of the (UUio et al. 2002) 
model: NFW Halo proflle with substructures (lower curves) and 
no-substructures (higher curves). The curves are the sensitivities 
considering the auto-correlation spectra information alone. 



I' indicates the multipoles bins. Let us also define /\(7*J''"" 



as the covariance matrix of the C^^g, 



covariance between C^^^j 



and C^Qgi, 



EGBl' 

i.e. AC^^", is the 
We divide the mul- 
tipole range 2 — 512 into 18 logarithmically spaced multipole 
bands. Adding to the previous signal, the contribution for 
a given DM model, specified through the annihilation chan- 
nel, the velocity averaged annihilation cross section and the 



DM mass. 



we further obtain C^^mA^'^J^'") ■ 



(.) which can 



be compared with the no-DM hypothesis as: 

x'^{{aAv),mx) = 

^i'ij,mn y-^EGBl' ~ ^DMl'J \^^ EGBV ) K'-^EGBV 



-CTmi') (11) 



where (AC;^"^™^ 



is the inverse of the covariance matrix. 
As a first step, however, we might want to consider a simple 
expression, for example considering as observables only the 
auto-correlation spectra C^egbi'- Unfortunately, the naive 
expression 
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turns out to be incorrect because the two spectra C^gbi' 
and C^J^n,! are correlated, i.e. AC!^;^7tI?, ^ 0. The correla- 



tion, indeed, is given by AC'^'aBV = 2/(2^ + l)(C|'£s^)^ i-e- 
by the cross correlation between the bands i and m, as intu- 
itively expected. The variance of Cegbi' can be also easily 
calculated and gives AC'^'aBv = 2/(2Z + l)iCiiaBi' + Cnf, 
where Cjv = 47r/A'evonts is the poisson shot noise. The cor- 
rect expression is then 
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To take into account multipole binning and partial sky cov- 
erage, in the above expressions we have to substitute 
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(16) 



where Al is the number of binned multipole in the band 1', 
and SNcounts/SQ is the density of the counts per steradian. 
Notice that in the above formulae no corrections for the 
angular resolution (the exp(Z^cri,) term) are required since 
we are referring to raw, PSF uncorrected spectra. Finally, 
we disregard the information at I < 10 to minimize the 
effects of the exposure. Let us also notice that, after all, the 
value of x^ calculated with Eg 1121 turns out to be generally 
similar to the one of Ea ll3l This behavior may be due to 
the different levels of noise in the various energy bands and 
in the autocorrelation spectra which makes the off-diagonal 
components in the covariance matrix (which do not have 
shot noise) generally subdominant. 

In case we want to exploit all the available informa- 
tion, i.e. including both the auto-correlations and the cross- 
correlations as observables, the full expression Eq llll for the 
X^ has to be employed. In this case the full covariance matrix 
becomes more complicated and we give the corresponding 
expressions in appendix B. 

Some further comments are in order regarding the cal- 
culation of x^ ■ We use always the same set of DM anisotropy 
maps to calculate the Cis for the cosmological DM case, de- 
spite the fact that they in principle depend on m^ (and the 
annihilation mode) through Eqs[6] and [T] Consequently, a 
different set of maps should be generated for each value of 
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m^ and for each channel. Fortunately, this dependence is 
weak above about 10 GeV where the DM horizon does not 
depend much on the DM annihilation energy spectrum so 
the effects on the sensitivity plots above 10 GeV are small, 
accordingly. To better justify this approximation we can also 
make a comparison with the case of DM annihilation from 
substructures in the Galactic halo where we use a very dif- 
ferent set of maps. The sensitivities generally change by a 
factor of a few in this case (see below). The slight depen- 
dence of the maps on DM model for the cosmological case 
is thus likely to be smaller than that. 

We can see the results in Fig. [T] where the sensitivities 
are reported for various annihilation channels. The abso- 
lute normalization of the DM spectra are obtained following 
the DM haloes clustering model of l|Ullio et al. 2002"|) as de- 
scribed in the introduction. Two versions of the model are 
shown: a conservative version with a NFW profile for the 
haloes and no substructures which gives a low DM normal- 
ization, and thus weaker constraints, and a more optimistic 
version where the DM signal is enhanced by the presence 
of substructures in the DM haloes. In principle, using the 
recent results of (|Zavala et al. 2009|) from the Millenium-II 
simulation, an order of magnitude enhancement with respect 
to the above "optimistic" case can be achieved with a corre- 
spondingly better sensitivity (see also (|Abdo et al. 2010a)) . 
in particular Fig.l). 

The various channels have generally similar behavior. 
The bb channel produces quite smooth limits. The fi'^ n~ 
and T^T~ sensitivity curves in contrast have a steeper slope 
and exhibit more structure. This is due to the fact that for 
these channels the gamma emission is concentrated at higher 
energies in a narrow peak near the energy corresponding 
to the DM mass. The sensitivity is thus more sensitive to 
the coarse binning in energy chosen for this analysis. The 
sensitivity in the t'^t~ channel and especially the fi^ fi~ 
channel at higher DM masses somewhat decreases due to the 
lower number of photons per annihilation (and thus lower 
statistics) available. 

Including the cross-correlation spectra information does 
not seem to improve the sensitivity. This probably has to 
be ascribed to the fact that the x^ for the auto-correlation 
case already includes some cross-correlation information in 
the off-diagonal terms of the covariance matrix. The level 
of anisotropies is also very low at low energy for this case 
(see Fig|3} and this somewhat reduces the importance of the 
information contained in the cross-correlation. The situation 
is different for the Galactic DM scenario as discussed below. 

We also checked the case in which only the highest en- 
ergy bands are employed for the analysis. We would expect 
an improvement in sensitivity given that very little DM sig- 
nal is expected in the 0.4-0.8 GeV band and thus scan- 
ning in this energy range introduces a statistical penalty 
factor. We find that the sensitivity indeed improves, but 
only marginally. On the other hand, the lowest energy band 
is useful because being mostly clean of DM it represents a 
good calibration for the background anisotropies. 

In Fig. |H]we show the case in which the DM anisotropy 
pattern does not vary with energy. This case is supposed 
to mimic the DM emission from unresolved substructures in 
the Milky Way. We thus simulate a single anisotropic map 
which we use for all energies. The angular spectrum for DM 
employed to simulate the map is normalized as in the case of 
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Figure 8. Same as previous figure but for the Galactic substruc- 
tures case. For a direct comparison with the cosmological case we 
use for the DM signal the same normalization(s) employed for the 
extra-galactic case. The solid curves are the sensitivities consid- 
ering the auto-correlation spectra information alone. The dotted 
ones are derived considering also the information contained in the 
cross-correlation spectra. 



substructures with a minimum mass of 10 Mq as derived in 
( [Siegal-Gaskins 2008[ l and it roughly corresponds to the ex- 
tragalactic anisotropies for an energy of about 300 GeV (see 
section 2). For a direct comparison with the cosmological 
case we use for the Galactic signal the same normalization 
employed for the extragalactic case, which is justified by the 
fact that more physically motivated models for the DM emis- 
sion in the MW have a large degree of freedom (|Ando 2009[) 
and can easily be tuned to match the extragalactic case. 
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The inclusion of the information from the cross- 
correlation spectra in this case, contrary to the cosmological 
one, gives a significant improvement, especially for higher 
DM masses. The reason for this difference is likely to be 
found in the very different set of anisotropy maps employed 
for the Galactic substructures case. In particular, the very 
high degree of anisotropy present also at low energy for this 
case enhances the importance of a direct measurement of 
the cross-correlation between low and high energy bands. 
In contrast, the very low DM anisotropy at low energy for 
the cosmological case (lower than the astrophysical EGB 
anisotropy) makes this information almost irrelevant. We 
also find in this case that excluding the lowest energy band 
does not improve the sensitivity significantly. 

The sensitivity increases substantially with respect 
to the cosmological case, especially for high DM masses. 
The slope of the sensitivity curves is indeed slightly less 
steep. Assuming higher levels of anisotropy for the Galactic 
case, which are possible, (|Ando 20091 |Siegal-Gaskins 2008| 
|Siegal-Gaskins and Pavlidou 2009| ) a correspondingly higher 
sensitivity could be achieved. However, the scaling is not lin- 
ear. It is worth noting for example that although at about 
10 GeV the Galactic anisotropics are about 2 orders of mag- 
nitudes larger than the extragalactic ones, the sensitivity is 
only a factor 2-3 better. The scaling however is better for 
m^ > 100 GeV where the difference in sensitivity is roughly 
a factor of 10. This non linear scaling is due to the lim- 
ited limited amount of detected photons (as it is the case at 
the edges of the sensitivity curve). Here, it is very difficult 
to distinguish between the two (galactic and extra-galactic) 
levels of anisotropics even if the latter is a factor 100 larger, 
since the statistical error bars are also very large. The much 
larger galactic anisotropy thus only modestly decrease the 
minimum amount of events required to detect the signal. 

Finally in Fig. [9] we show the impact on the DM sensi- 
tivity of changing the underlying astrophysical EGB model 
to the Poisson dominated one. For illustration, only the ex- 
tragalactic T^T^ case with NEW profile and substructures 
is shown. The interesting results (perhaps somewhat counter 
intuitive) is that even with this increased level of background 
in the anisotropy spectrum, still the sensitivity to the DM 
component remains roughly unchanged. The reason is that a 
Poisson noise at the level of 5.0e-5 is anyway clearly detected 
with a statistics of 5 years, so the DM component can still 
be easily detected as a modulation (as usual, in multipole 
and in energy) of this dominant Poisson noise anisotropy 
spectrum. The underlying assumption here, however, is that 
this new background is fully under control and well charac- 
terized. The extent to which this will be true is however 
not clear and a better understanding may require dedicated 
simulations and analyses which are beyond the scope of this 
paper. 



8 SUMMARY AND CONCLUSIONS 

In the present analysis we have calculated the sensitivity 
to Dark Matter annihilation exploiting the characteristic 
anisotropy properties that DM would imprint in the Extra- 
Galactic Gamma Background. We performed a realistic sim- 
ulation of the Fermi-LAT including the effect of the Point 
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Figure 9. DM sensitivity for annihilation into into 66, extragalac- 
tic anisotropics and NFW Halo profile with substructures, for the 
two cases of EGB background anisotropics from the correlation 
term only (same case of Fig|8l solid contours) and from a pure 
Poisson term only (dashed contours). 



Spread Function, effective area and of the anisotropic sky 
exposure. 

More importantly, we investigated the impact of Galac- 
tic foregrounds on the DM anisotropy sensitivity using a 
GALPROP model to produce the sky maps of the Galactic 
gamma emission. The contribution of Galactic foregrounds 
is generally non negligible even above Galactic latitudes of 
about 60°. Foreground cleaning will thus be an important 
(and likely challenging) issue in the analysis of the real data. 
For the present work we do not attempt foreground cleaning 
(except for masking the Galactic plane and point sources) 
and we only require that the anisotropy from DM has to ex- 
ceed the level of foregrounds anisotropy to be detectable. 
Further, we only analyze a high Galactic latitude region 
where the model dependence of the foreground is expected 
to be moderate. In this respect our analysis can be regarded 
as a conservative one: In reality, if foreground cleaning can 
be effectively implemented, a sensitivity to a lower level of 
anisotropy can be achievable. Also, a larger fraction of the 
sky can be employed than the small conservative one used in 
the present analysis. We leave a more careful investigation 
of the issue for future work. 

Despite the difficulties involved in the foreground anal- 
ysis, we find that the prospects for detecting DM anisotropy 
signatures are promising. For some WIMP masses, sen- 
sitivity to {(7A«) ~ 3 X 10~'^®cm^s~^(i.e. the thermal 
relic value) can be achieved using the anisotropy infor- 
mation alone. Interestingly, these constraints are compara- 
ble to what is achievable using the energy spectrum alone 
(|Abdo et al. 2010a)) . The anisotropy information is thus 
complementary to the energy spectrum and provides a cross- 
check and an increase in the overall sensitivity to the DM 
signal. 

We also investigated the recently proposed scenario of 
a population of DM substructures in the Milky Way halo as 
predicted by N-body simulations. These subhaloes give rise 
to an annihilation signal which is almost uniform in the sky 
and which can be comparable in intensity to the observed 
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isotropic gamma background. For this case, the anisotropy 
signal is predicted to be even higher than for the extragalac- 
tic case due to the huge dumpiness of the DM distribution 
that correspondingly gives an enhanced sensitivity. We find, 
indeed, that when increasing the anisotropy the sensitivity 
to the DM signal increases as well with an improvement of 
a factor of ~10 for m^ > 100 GeV. 

In summary, simulations indicate that ,with the sensi- 
tivity of the Fermi-LAT, detailed anisotropics studies of the 
gamma-ray sky can be performed. The next obvious step 
will be the analysis of the actual data. Anisotropy analyses 
will soon provide a complementary source of information 
which will give further insights into the DM problem and 
the properties of high-energy astrophysical sources. 
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where < ... > indicates the statistical ensemble average. The 
quantity 



Ci = 



v^ air, 



^^21 + 1 



(A2) 



is an unbiased estimator of the true power spectrum Ci, i.e. 
<C, >=Ci. 

In the case of two maps, besides the two autocorrelation 
power spectra C} and Cf , further information is contained in 
the cross- correlation power spectrum which can be estimated 
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where the equivalence of the two expressions comes from 
the equality: a\^ = (— l)'"a;1mi i = 1,2. The three spectra 
are of course not independent but are part of a covariance 
matrix: 



cl Cl 



(A4) 



which has to be talcen into account when evaluating the 
errors for the forecast. This is described in more details in 
the next section. 

In practice, to calculate the power spectra, the 
harmonic transform aj^ is performed with HEALpix 
(|Gorski et al. 20051 after binning the events into maps of the 
HEALpix format itself. Due to the finite number A'' of col- 
lected events beside the intrinsic anisotropies the data also 
contain noise (shot noise or white noise) which appears in 
the power spectrum as a constant whose value is 4n/N. The 
noise sets the limit above which the intrinsic anisotropies can 
be detected. It is standard to remove the shot noise from the 
spectrum showing only the intrinsic power spectrum. It is 
also customary to show not the power spectrum Ci but the 
quantity l{l + l)Ci/2n. 



APPENDIX B: 
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COVARIANCE MATRIX 



In the general case in which we include as observables the 
cross-correlation power spectra the covariance matrix can 
be calculated by propagation of the errors from the expres- 
sions |X2] and [XS] and then taking the ensemble average. The 
variance AC^^^j, of the cross-correlation spectrum C^^„g^, 
{i 7^ j) can be written as 



21 
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APPENDIX A: 
SPECTRUM 



ANGULAR POWER 



We summarize in this appendix the conventions and nota- 
tion used for the anisotropy analysis. The anisotropy of a sig- 
nal on the sky can be quantified in terms of the power spec- 
trum i.e. the spherical harmonic transform on the sphere. 
]VIore specifically, a map on the sphere f{Q) can be decom- 
posed in spherical harmonics as f{Cl) = E;m,a;m,y^(f2) and 
the angular power spectrum can then be defined as 



We then have the case of covariance AC^J^, between 
an auto-correlation spectrum C^^Qgy and cross-correlation 
spectrum C^q^i, . For this case we have 
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For the case of covariance between two cross-correlation 
spectra we have to distinguish two cases. In the case of co- 
variance AC'^Qgi, between C^^^gi, and C^'Jjgj, where there 
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is an energy band (an index) in common we have: 

^^EGBV —^ °'^EGBl'"'^EGBl' > - 

27^ {(C'iGBV + C^d^C^GBv) + {C%^bi:){C'egbv)) (B3) 

If instead there is no energy band (no index) in common 
^^^egbv finally we have 

^'■^EGBl' -< Oi^EGBl'°^EGBl' > — 

^ {{CTgbv){C'e"gbv) + {C^EGBv){C%^GBi')) (B4) 

Expression (|14|) and (IBlf) can also be rewritten in the 



more familiar form 



5Ci _ / 2(1 + Civ/Ci) 



Ci V (2' + l)AZ/.ky ' 

for the auto-correlation spectra and by 



(B5) 




(B6) 



for the cross-correlation spectra, where again 

Civ. = (^^%p) /sky (B7) 

and AZ is the number of binned multipole in the band V , 
and (5 Accounts /^^ is the density of the counts per steradian 
in the map i. 



